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Abstract 

We explore the phenomenology of Kaluza-Klein (KK) dark matter in very general models with 
universal extra dimensions (UEDs), emphasizing the complementarity between high-energy collid- 
ers and dark matter direct detection experiments. In models with relatively small mass splittings 
between the dark matter candidate and the rest of the (colored) spectrum, the collider sensitivity 
is diminished, but direct detection rates are enhanced. UEDs provide a natural framework for such 
mass degeneracies. We consider both 5-dimensional and 6-dimensional non-minimal UED models, 
and discuss the detection prospects for various KK dark matter candidates: the KK photon 71, 
the KK Z-boson Z\, the KK Higgs boson H\ and the spinless KK photon jh- We combine collider 
limits such as electroweak precision data and expected LHC reach, with cosmological constraints 
from WMAP, and the sensitivity of current or planned direct detection experiments. Allowing 
for general mass splittings, we show that neither colliders, nor direct detection experiments by 
themselves can explore all of the relevant KK dark matter parameter space. Nevertheless, they 
probe different parameter space regions, and the combination of the two types of constraints can 
be quite powerful. For example, in the case of 71 in 5D UEDs the relevant parameter space will 
be almost completely covered by the combined LHC and direct detection sensitivities expected in 
the near future. 

PACS numbers: 95.35.+d,11.10.Kk,12.60.-i,95.30.Cq,95.30.-k,14.80.Ly 
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I. INTRODUCTION 



The Standard Model (SM) has been extremely successful in explaining all available ex- 
perimental data in particle physics. However, there are several unsettling features of the SM, 
which have motivated a substantial research effort on physics beyond the Standard Model 
(BSM). The two issues continuously attracting the most attention are the hierarchy problem 
and the dark matter problem. The anticipated discovery of the Higgs boson of the SM at 
the Large Hadron Collider (LHC) at CERN would pose a challenging theoretical question: 
what is the next fundamental energy scale? If it is as high as the Planck scale, then what 
stabilizes the hierarchy between the Planck and electroweak scales? Or, if it is much lower 
than the Planck scale, what is the physics associated with it? The second issue is related 
to the now established existence of a dark matter (DM) component of the universe. Since 
the SM does not accommodate a suitable DM particle candidate, the dark matter problem 
is the most pressing phenomenological evidence for physics BSM 

A. The Dark Matter Problem and Physics Beyond the Standard Model 

There are different avenues one could follow in extending the SM and addressing the 
dark matter problem. The common theme among them is the introduction of new particles, 
one of which is neutral and serves as the dark matter candidate; and a new symmetry, a 
remnant of which survives in the low energy effective theory and ensures that the lifetime 
of the DM particle is sufficiently long (at the minimum, longer than the age of the uni- 
verse). In principle, simply postulating a new stable and neutral particle would be rather 
ad hoc and unsatisfactory without further corroborating evidence. Fortunately, the DM 
candidates in most BSM models typically have some kind of non-gravitational interactions, 
which are sufficient to keep them in thermal equilibrium in the early universe. Thus, their 
relic abundance can in fact be straightforwardly calculated in any given model (for details, 
see Section Til Bl below) . The generic result of this computation is that a weakly interacting 
massive particle (WIMP) with a mass near or below the TeV scale has a relic density in 
the right ballpark, and is a suitable candidate for dark matter. By now there are many 
examples of WIMPs in BSMs, perhaps the most popular being the lightest superpartner 
(LSP) in supersymmetry (SUSY) with R-parity conservation [2j, the lightest Kaluza-Klein 
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partner (LKP) in Universal Extra Dimensions |3|, the lightest T-parity odd particle in Little 

Tin n n 

Higgs models |4|, |5|, the lightest U-parity odd particle in U(l) '-extended models |6|, |7|, etc. 

The most exciting aspect of the WIMP DM hypothesis is that it is testable by experiment. 
Indeed, WIMPs near the TeV scale can be easily within reach of both high-energy colliders 
and dark matter detection experiments. Furthermore, the size of the corresponding DM 
signals can be readily calculated within any given BSM, providing some rough expectations 
for discovery in each case. In principle, the signals depend on a typically a large number of 
model parameters. However, speaking in a broader sense, the WIMP DM phenomenology 
mostly depends on the answers to the following two questions: 

• Ql: What is the identity of the DM particle candidate? 

• Q2: What is the size of the mass splitting between the DM particle and the rest of 
the (relevant) spectrum? 

In the following two subsections we shall discuss each one of these questions and thus moti- 
vate our setup and methodology. 

B. The Nature of the Dark Matter Particle 

Within any given BSM, there are typically several potential dark matter candidates (i.e. 
neutral and stable particles) present in the spectrum. The answer to the first question (Ql) 
therefore selects one of them as the "true" dark matter. For example, in SUSY, the dark 
matter particle could be either a fermion (e.g. gravitino or the lightest neutralino) or a boson 
(the lightest sneutrino). In turn, the lightest neutralino could be the superpartner of a gauge 
boson (e.g. a Bino, a Wino, possibly a Z'-ino), the superpartner of a Higgs boson (e.g. a 
Higgsino or a singlino) or some admixture of these j^j]. Similarly, the lightest sneutrino could 
carry anyone of the three lepton flavors, and in addition, could be left-handed js], right- 



handed [9J, or some mixture of both 10j. Since all of these particles have rather different 
properties, it is clear that it is impossible to make any generic predictions about SUSY dark 
matter without specifying the exact nature of the LSP, i.e. providing the answer to Ql 
above. 

On the positive side, the answer to Ql goes a long way towards the determination of 
the size of the expected dark matter signals. Once the identity of the dark matter particle 



is specified, its couplings are fixed and can be used in the calculation of both direct and 
indirect detection rates. What is even better, the answer to question Ql can be provided 
in a rather model-independent way, without reference to the exact specifics of the model, 
such as the physics of the ultraviolet completion, Renormalization Group Equation (RGE) 
evolution down from high scales, etc. 

In this paper, we shall explore the dark matter phenomenology of general models with 



flat universal extra dimensions llj , where the usual Standard Model structure is embedded 
in 5 or 6 space-time dimensions. We shall assume the same gauge symmetry and particle 
content as in the SM. Similar to the SUSY case just discussed, the models contain several 
possible dark matter candidates (electrically- neutral particles which are stable due to KK 
parity conservation) 1 . In five dimensional models with minimal particle content, they are: 
the KK graviton (G\), the KK neutrino (i^i), the KK photon (71), the KK Z-boson (Z\) and 
the KK Higgs boson (Hi). Six dimensional UED models present additional possibilities: the 
spinless KK photon (7^) and the spinless KK Z-boson (Zh), which are linear combinations 
of the gauge boson polarizations along the two extra dimensions. Just like the case of 
SUSY, which of these particles is the lightest and thus the dark matter candidate, depends 
on the model-building details. The issue is even more subtle than in SUSY, since all of 
these KK particles have tree-level masses of the same order, proportional to the inverse 
radius R~ l of the extra dimension. This mass degeneracy is lifted by two main sources: 
radiative corrections due to renormalization and nonuniversality in the boundary conditions 
at the cut-off scale. The former effect is in principle computable within any given model 

nnn 

PJ, H3|, Hi|, while the latter is a priori unknown, as its origin lies in the ultraviolet physics 



above the cut-off scale 



0, Q, 



16 



171 ]. A common assumption throughout the existing 
literature on UED is to ignore any boundary terms at the cut-off scale. The resulting model 
has been dubbed "Minimal UED" and is known to accommodate only ji and Gi LKP in 
five dimensions lj], Ji| and 7# in six dimensions [lj], [^J. However, given our complete 
ignorance of the physics at and above the cut-off scale, the other possibilities for the nature 
of the KK dark matter particle should be given serious consideration as well. 

One of the goals of this paper is to start filling these gaps in the literature, by exploring 
the phenomenology of the alternative dark matter candidates in UED. Of course, not all of 



1 For further details on UED models, see Section Hi At 
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them are on an equal footing. For example, the KK graviton G\ interacts with the Standard 
Model particles too weakly to be relevant for direct detection searches. The KK neutrino V\ 
is already ruled out due to its large elastic scattering cross section [21] . We shall therefore 
concentrate on the remaining two possibilities in 5D UEDs: the KK Z-boson (Z\) and the 
KK Higgs boson (Hi). We shall also review and update the previously published results on 
71 and 7#, so that our work would provide a concrete and complete reference on KK dark 
matter. 



C. The Effect of a Mass Degeneracy on Dark Matter Signals 

The second important issue for dark matter phenomenology is the answer to question 
Q2, namely, what is the mass splitting between the dark matter particle and the rest of 
the spectrum. Of course, it is in principle possible to have the dark matter particle as the 
only new particle in the model, in which case Q2 does not apply, and the predictions for 
the dark matter signals are quite robust, once Ql is addressed. However, realistic models 
typically contain a multitude of new particles, in addition to the dark matter candidate. 
Their proximity (in mass) to the dark matter particle therefore becomes an important issue, 
at least in three, very different aspects. 

The first is related to the predicted dark matter relic abundance. A close mass degeneracy 



can increase the importance of coannihilation processes at freeze-out 22| . and the results 
for the relic density are now sensitive not only to the properties of the dark matter particle 
itself, but also to the properties of the coannihilating particles. The size of the coannihilation 
effect depends on the particular scenario, and there are several known cases in which it can 
be significant, e.g. Bino-like neutralinos in supersymmetry. The calculation of the relic 
density in the presence of coannihilations is a bit more involved (due to the larger number of 
processes which need to be considered), but nevertheless pretty straightforward. For UED 
models, where mass degeneracies are generically expected, the complete set of coannihilation 
orocesses which are relevant for the 71 and Z\ LKP case in 5D UED have been calculated 



21 



23, 



24j . We shall make use of them in our analysis below in Section Hi B[ After reviewing 
the case of 71 LKP, which has been previously discussed in the context of minimal UED, we 
shall also consider Z\ LKP and illustrate the effects of coannihilations with KK quarks on 
its relic abundance. Since a calculation of coannihilations in 6D UED models is still lacking, 
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there we shall consider only one specific example in detail - the previously discussed case of 



251 ] . The corresponding results for the direct detection rates of Zh can be obtained by 



a simple scaling of the gauge couplings. 

A small mass splitting also has a large impact on the expected direct detection signals, 
whenever the particle degenerate with the LKP can be exchanged in an s-channel. This 
situation may in principle arise in supersymmetry, if the squarks are very light, but this 
would be viewed by most people as a fortuitous accident. On the other hand, such a 
degeneracy occurs much more naturally in UED, where the masses of the KK quarks and 
the LKP necessarily have a common origin (the scale of the extra dimension). The mass 
degeneracy may lead to a substantial enhancement of the LKP elastic scattering rate [26] . In 
Section III CI we first review the calculation of the spin-independent and the spin-dependent 
elastic scattering cross sections for the 71 LKP case 26j. Then we also consider the case of 



Z\, Hi and 7^ LKP, paying special attention to the enhancement of the cross sections in 
the limit of small mass splittings. 

Finally, the mass splitting between the dark matter candidate and the rest of the new 
physics spectrum is an important parameter for collider searches as well. The discovery 
reach for new physics at colliders is greatly diminished if the mass splittings are small. 
This is because the observable energy in the detector would then be rather small as well, in 
spite of the large amount of energy present in the events. Correspondingly, the measured 
missing energy (and any related variable such as Ht) is also rather small, which makes it 
more difficult to extract the new physics signal from the SM backgrounds. Fortunately, as 
mentioned above, this is precisely the case when direct detection is more promising. In Sec II VI 
we shall explore this complementarity for various KK DM scenarios, focusing on KK gauge 
boson dark matter. From the previous discussion it should be clear that having specified 
the nature of the DM particle, the two most relevant parameters are the DM particle mass 
m x and the mass splitting with the nearest heavier colored particles. In Sec |IV] we shall 
utilize this two-dimensional parameter space, and contrast constraints from different sources: 
colliders, cosmological observations, and current or planned direct detection experiments 
(the latter are first extensively reviewed in Sec. IHIj) . As expected, we find that colliders and 
dark matter searches are highly complementary, while the WMAP constraint is orthogonal 
to them but is somewhat more model-dependent. Section |V] is reserved for a summary and 
conclusions. In the appendix we write out some technical details of our analysis. 
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II. UNIVERSAL EXTRA DIMENSIONS AND KALUZA-KLEIN DARK MATTER 



A. Review on Universal Extra Dimensions 

Models with universal extra dimensions place all Standard Model particles in the bulk of 
one or more compactified flat extra dimensions. In the simplest and most popular version, 
there is a single extra dimension compactified on an interval, Si/Z 2 . In UED, each SM 
particle has a whole tower of KK modes. The individual modes are labelled by an integer 
n, called KK number, which is nothing but the number of quantum units of momentum 
which the SM particle carries along the extra dimension. A peculiar feature of UED is 
the conservation of Kaluza-Klein number at tree level, which is a simple consequence of 
momentum conservation along the extra dimension. 

However, the fixed points in orbifold compactifications break translation invariance along 
the extra dimension. As a result, KK number is broken by bulk and brane radiative ef- 
fects 0, Q, Q down to a discrete conserved quantity, the so called KK parity (— l) n . 
The geometrical origin of KK parity in the simplest (Si/Z 2 ) case is the invariance under 
reflections with respect to the center of the interval. Since KK parity is conserved, the 



lightest KK parity odd particle is a suitable WIMP candidate [14], [21|, [26|, [27]. KK parity 
also ensures that the KK-parity odd KK partners (e.g. those at level one) are always pair- 
produced in collider experiments. This is reminiscent of the case of supersymmetry models 
with conserved .R-parity. Therefore, the limits on UED KK modes from collider searches 
are relatively weak and are rather similar to the limits on superpartners. KK-parity is also 
responsible for weakening the potential indirect limits on UED models from low-energy pre- 
cision data. Just like SUSY models with R-parity, the virtual effects from new physics only 



28, 29, 3 



appear at the loop level and are loop suppressed 

Since all KK modes carry momentum along the extra dimension, at tree-level their masses 
receive a dominant contribution -| ~ n TeV, and a subdominant contribution from the 
corresponding SM particle mass. All KK modes at a given KK level n are therefore quite 
degenerate. The KK modes of the lightest SM particles (photons, leptons, light quarks) 
even appear to be absolutely stable at tree level. However, this conclusion is invalidated 
after accounting for the radiative corrections to the KK masses. The latter are proportional 
to |r and are sufficient to lift the degeneracy between the lightest KK modes, leaving only 



S 



one of them (the true LKP) as absolutely stable jl^j ]. 

The nature of the LKP, on the other hand, is more model-dependent. In the minimal 
5D UED model, where the boundary terms at the cut-off scale are ignored, the lightest 
KK particle is typically the n = 1 mode Bi of the hypercharge gauge boson [141 ] . Since 
the Weinberg angle for the level one neutral gauge bosons is rather small, B\ is essentially 
also a mass eigenstate, the KK "photon", and we shall therefore denote it as 71. The 
KK photon 71 is an attractive dark matter candidate 2l|, |26j, whose relic abundance is 
consistent with the observed dark matter density for a mass range between 500 GeV and 
about 1.5 TeV, as shown by detailed computations including coannihilations [23, [24j and 



level-2 resonances 



31 



32l . |33j . Direct detection of this KK dark matter may be within reach 



of the next generation experiments 



26 



35 



361 ] . Indirect detection of K 



^ dark matter also 
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40 



41 



42 



43 



44] 



has better prospects than the case of neutralinos in SUSY 

In UED the bulk interactions of the KK modes readily follow from the Standard Model 
Lagrangian and contain no unknown parameters other than the mass of the Standard 
Model Higgs boson. In contrast, the boundary interactions, which are localized on the 
orbifold fixed points, are in principle arbitrary, and thus correspond to new free parameters in 



the theory. They are in fact renormalized by bulk interactions, and are scale dependent [12]. 
Therefore, we need an ansatz for their values at a particular scale. Virtually all existing 
studies of UED have been done within the framework of minimal UED (MUED), in which 
the boundary terms are assumed to vanish at the cut-off scale A, and are subsequently 
generated through RGE evolution to lower scales (see for 5D and HQ for 6D). 

In the minimal UED model therefore there are only two input parameters: the size of the 
extra dimension R, and the cut-off scale A. Of course, there are no compelling reasons 
for assuming vanishing boundary terms: the UED model should be treated only as an 
effective theory which is valid up to the high scale A, where it is matched to some more 
fundamental theory, which is generically expected to induce nonzero boundary terms at the 
matching scale A. As already mentioned in the introduction, nonvanishing boundary terms 
may change both the nature of the LKP, as well as the size of the KK mass splittings. The 
resulting phenomenology may be very different from the minimal case. This is why in this 
paper we shall allow for more general scenarios with Z\ and Hi LKP 2 . In each case, we shall 



Other dark matter candidates, such as the level-1 KK mode of the graviton or of a right-handed neutrino, 
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take the LKP mass tulkp and the LKP - KK quark mass splitting 



_ m qi - m LKP 

91 ~~ ™ ' V ' 

as free parameters. We remind the reader that after compactification, the low energy effective 
theory contains two massive (Dirac) KK fermions for each (Dirac) fermion in the Standard 
Model. The KK fermions are properly referred to as S'f/(2)ty-doublet KK fermions or 
S , f/(2)iy-singlet KK fermions. However, in the literature they are sometimes called "left 
handed" and "right handed" , referring to the chirality of the corresponding Standard Model 
fermion at the zero level of the KK tower. This nomenclature may lead to some confusion, 
since all KK fermions are Dirac and have both chiralities. In our study, we shall treat the 
SU (2)n/-doublet KK quarks (often denoted by Q\) and the SU (2) ^-singlet KK quarks (often 
denoted by qi) equally, thus avoiding the need for two separate mass splitting parameters 
(for example, a separate Aq 1 and A 9l ). The generalization to the case of different KK quark 
masses is rather straightforward. 

We shall also explore cases with more than one universal extra dimension. Theories 
with two universal extra dimensions also contain a KK parity. Under the simplest com- 
pactification which leads to chiral zero- mode fermions (a "chiral" square with adjacent sides 



identified 



51 



521]). the KK parity transformations are reflections with respect to the cen- 
ter of the square. Momentum along the two compact dimensions is quantized so that any 
6-dimensional field propagating on the square appears as a set of 4-dimensional particles 
labelled by two positive integers, (n,m). These particles are odd under KK parity when 
n + m is odd and are even otherwise. In any process, odd particles may be produced or 
annihilated only in pairs. The lightest odd particle, which is one of the (1,0) states, is 
thus stable. Gauge bosons propagating in six dimensions may be polarized along the two 
extra dimensions. As a result, for each spin-1 KK particle associated with a gauge bo- 
son, there are two spin-0 KK fields transforming in the adjoint representation of the gauge 
group. One linear combination becomes the longitudinal degree of freedom of the spin-1 
KK particle, while the other linear combination remains as a physical spin-0 particle, called 



are also viable for certain ranges of parameters in models with one universal extra dimension [18l . |46 



47 



10 



the spinless adjoint 3 . In the minimal model with vanishing boundary terms, the radiative 



corrections 



3. 



201 ]. are such that the lightest (1,0) particle on the chiral square [5l|, |52J is 



always a linear combination of the electrically-neutral spinless adjoints of the electroweak 
gauge group. Due to the small mixing angle, this linear combination is essentially a photon 
polarized along the extra dimensions. Similar to its_5D cousin 7-1 , the sp inless photon 7^ 
in 6D UED is also a viable dark matter candidate 



25] . See Refs. 



53 



54. 



55| for KK dark 



matter candidates in UED models with an extended gauge symmetry. 



B. Relic Density Calculation with Coannihilations 

We briefly review the calculation of the relic density including coannihilation processes. 
When the relic particle x is nearly degenerate with other particles in the spectrum, its relic 
abundance is determined not only by its own self-annihilation cross section, but also by 
annihilation processes involving the heavier particles. The generalization of the relic density 



calculation including this "coannihilation" case is straightforward |21l . |22j . Assume that the 
particles x% are labelled according to their masses, so that m ; < nij when % < j. The number 
densities rii of the various species Xi o bey a set of Boltzmann equations. It can be shown 



that under reasonable assumptions [22|, the ultimate relic density n x of the lightest species 
Xi (after all heavier particles Xi have decayed into it) obeys the following simple Boltzmann 
equation 

dn 



-jf = ~3Hn x - (<T ef fv)(n 2 x - n 2 eq ) , (2) 

where H is the Hubble parameter, v is the relative velocity between the two incoming 
particles, n eq is the equilibrium number density and 

JV 

a eff (x) = J] C T 4J ^(l + A J ) 3 / 2 (l + A J ) 3 / 2 exp(-x(A J + A J )) , (3) 

ij g eff 
N 

9eff(x) = ^&(1 + A i ) 3 / 2 exp(-a;A i ) , (4) 

i=l 

m i -m 1 mi 
mi 1 



3 



In contrast to the 6D case, KK particles in 5D UED arc labelled by only one integer and spinless adjoints 
do not exist since there is only one gauge boson degree of freedom polarized along the extra dimension. 
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Here cry = a(xiXj ~^ SM) are the various pair annihilation cross sections into final states 
with SM particles, gi is the number of internal degrees of freedom of particle \% an d n x = 
J2iLi n i * s the density of Xi we want to calculate. 



By solving the Boltzmann equation analytically with appropriate approximations 
the abundance of the lightest species Xi is given by 



21 



22|, 



2 1.04 x 10 9 GeV" 1 x F 1 



X M « y/gJ^jla + Mb/xp' 

where the Planck mass scale is Mpi = 1.22 x 10 19 GeV and g* is the total number of effectively 
massless degrees of freedom at temperature T: 

g.(T)= Yl 9i + l £ ( ? ) 

i=bosons i=fermions 

The functions J a and /b are defined as 

P CO 

I a = %f a e jf(x)x~ 2 dx , (8) 

J Xp 

POO 

h = 2x 2 F / b e ff(x)x~ 3 dx . (9) 

J Xp 



The freeze-out temperature, x F , is found iteratively from 

Xb _ ln ^ c ( c | o) f^ 9eff(xF)m 1 Mp l (a eff (x F ) + 6b ef f(x F )/x F ) \ 
V V 8 2tt 3 y/g*{x F )x F ) 

where the constant c is determined empirically by comparing to numerical solutions of the 
Boltzmann equation and here we take usual. a e ff and b ef f are the first two terms 

in the velocity expansion of cr e // 

°Vf/(z) ^ = a e// 0) + &e//(z) v 2 + 0(> 4 ) . (11) 



12 



Comparing Eqns. (j3J) and (jTTI) . one gets 



N 



aeff(x) = ^a,^(l + A l ) 3/2 (l + A i ) 3/2 exp(-a:(A J + A J )) 



N 



9eff 



5>,^(1 + A,) 3 / 2 (l + A J ) 3/2 exp(-x(A, + A,)) 



(12) 
(13) 



where and fey are obtained from a^v = + bijV 2 + 0{v 4 ) and t> is the relative velocity 
between the two annihilating particles in the initial state. Considering relativistic corrections 
56J to the above treatment results in an additional subleading term which can be accounted 
for by the simple replacement 

b^b-^a, (14) 



in the above formulas. For our calculation of the relic density, we use the cross sections 
given in Refs. 



21 



23 



24]. 



As explained earlier, the assumptions behind MU 
nonvanishing boundary terms at the scale A 



15 



16. 



17 



D can be easily relaxed by allowing 



45( | . This would modify the KK spec- 



trum and correspondingly change the MUED predictions for the KK relic density. Within 
the modified KK spectrum, any neutral KK particle could be a dark matter candidate. As 
an illustration here we shall consider the case of 71 and Z\ LKP 4 , for which the results for 
the relevant coannihilation processes are available in the literature 



23, 



2j. In Fig. [H we 



show the relic densities of 71 and Z\ as a function of the corresponding LKP mass (m 7l or 
Tnz 1 ) in 5D UED. We include coannihilation effects with all n = 1 KK particles with prop- 
erly defined masses. The (black) solid lines show the LKP relic density for several choices 
of the mass splitting (0Q) between the LKP and the KK quarks. We assume that singlet and 
doublet KK quarks are degenerate (i.e., Aq 1 = A qi ). The green horizontal band denotes the 
preferred 2<r-WMAP region for the relic density 0.1037 < QcDAih 2 < 0.1161 [581 ] . The cyan 



vertical band delineates values of rriLKP disfavored by precision data 



59l . |60j 5 . In each case 



of Fig. [T^,, we use the MUED spectrum to fix the masses of the remaining particles, and 



4 The Z\ is also a good dark matter candidate in warped extra dimensions with KK parity 57 1. 

5 While there have been no studies of non-minimal UED models with Z\ LKP, we anticipate that the 
precision bounds in that case will be similar to those in MUED, therefore we display the same indirect 
constraint in Fig. [TJ). 
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FIG. 1: Relic density of the LKP ((a) 71 and (b) Z\) as a function of the LKP mass. The (black) 
solid lines show the LKP relic density for several choices of the mass splitting |7]) between the LKP 
and the KK quarks. We assume that singlet and doublet KK quarks are degenerate. The green 
horizontal band denotes the preferred WMAP region for the relic density 0.1037 < ClcDAih 2 < 
0.1161. The cyan vertical band delineates values of tulkp disfavored by precision data, (a) We 
vary the qi mass by hand, keeping the masses of the remaining particles fixed as in MUED. The 
solid lines from top to bottom correspond to A qi = 0.5,0.1,0.05,0.02,0.01. The (red) dotted line 
is the result from the full calculation in MUED, including all coannihilation processes, with the 
proper MUED choice for all masses, (b) We assume Z\ and W^ 1 are degenerate, the gluon is 
heavier than Z\ by 20%, while all other KK particles are heavier than Z\ by 10%. The solid lines 
from top to bottom correspond to A qi = 0.01,0.02,0.05,0.1,0.5. 



then vary the (common) KK-quark mass m qi by hand. The solid lines from top to bottom 
correspond to A gi = 0.5, 0.1, 0.05, 0.02, 0.01. The (red) dotted line is the result from the full 
calculation in MUED, including all coannihilation processes, with the proper MUED choice 
for all masses. In Fig. [lb we assumed Z\ and W x are degenerate, the gluon is heavier than 
Zi by 20%, while all other KK particles are heavier than Z\ by 10%. The solid lines from top 
to bottom correspond to A qi = 0.01,0.02,0.05,0.1,0.5. Some individual quantities entering 
the relic density calculation for 71 (Z\) LKP are shown in Fig. [2] (Fig. [3]). 

We see that coannihilations in the case of 71 LKP decrease the prediction for Qh 2 and 
therefore increase the range of preferred m 7l values. For A 9l on the order of a few percent, 
the desired range of m 71 is pushed beyond 1 TeV. This poses a challenge for any collider 
searches for UED, since the KK production cross sections at the LHC become kinematically 
suppressed for heavier KK modes. What is even worse, the small mass splitting A 91 degrades 
the quality of the discovery signatures, e.g. the cascade decays of the KK quarks would yield 
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FIG. 2: Plots of various quantities entering the 71 LKP relic density computation, as a function 
of the mass splitting A qi between the LKP and the KK quarks: (a) relic density, (b) a e ff(xp), (c) 
g e ff(xp), and (d) xp. In all four panels, the KK quark masses are varied by hand according to 
the value of A qi , while the masses of the 71 and the remaining KK modes are held fixed at their 
nominal values predicted in MUED for R^ 1 = 1 TeV. 



only (rather soft) jets and no leptons. 

On the other hand, Fig. [Ib reveals that coannihilations with KK quarks have the opposite 
effect in the case of Z\ LKP 6 . This time the effect of coannihilations is to increase the 
prediction for Qh 2 and thus lower the preferred range of values for mz x . The lesson from 
Figs. [T^l and [T|d is that while coannihilations can be quite important, the sign of the effect 
cannot be easily predicted, since, as will be illustrated in Figs. [2] and [31 it depends on the 
detailed balance of several numerical factors entering the computation. We shall discuss 



6 A similar behavior exists in the case of 71 LKP when coannihilations are caused by the SU (2)-singlet KK 
leptons £ri 
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FIG. 3: The same as Fig. \2\ but for Z\ LKP. Here the Z\ and Wi masses are taken to be 2 TeV, 
the KK quark masses are varied in accordance with A qi , while the remaining spectrum is fixed as 
in Fig. [I]b. 



these in some detail in the remainder of this subsection. Readers who are not interested in 
these numerical details, are invited to jump to Section Hi CI 

In Fig. [2^ (Fig. [3^) we plot the relic density of the 71 [Z\] LKP, as a function of the mass 
splitting A qi between the KK quarks and the corresponding LKP. The rest of the spectrum is 
held fixed as explained in the figure captions. Figs.[2h and[3h demonstrate the importance of 
coannihilations at small mass splittings. For A qi larger than about 10 — 20%, coannihilations 
are turned off, but for KK quarks within 10% of the LKP mass, the coannihilation effect is 
significant. For 7x LKP, it lowers the prediction for the relic density Qh 2 , while in the case 
of Zi LKP Qh 2 is enhanced. In order to understand this different behavior, it is sufficient 
to investigate the coannihilation effect on the effective cross section, and in particular the 
dominant term a e //> which is plotted in Figs. [2b and[3b. As can be seen from eq. (TI2"|) . every 
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term contributing to a e // is a ratio between two quantities, each of which has a nontrivial 
A c/1 dependence. The denominator is common to all terms and is nothing but the effective 
number of heavy particle degrees of freedom g e ff defined in eq. We show the g e ff 

dependence on A qi in Figs. [2b and [3b. As expected, g e ff increases significantly after the 
turn-on of coannihilations (below A qi ~ 0.1), due to the large multiplicity of KK quark 
states. At the same time, the numerator of each term contributing to the a e ff sum ( Fl2|) is 
simply the Boltzmann suppressed annihilation cross section, which also increases with the 
onset of coannihilations (at small mass splittings A qi ). The net effect on a e ff is determined 
by which of these two quantities increases faster at small A qi , relative to the nominal case 
without coannihilations. In the case of 71 LKP, the self-annihilation cross sections are rather 
weak, due to the smallness of the hypercharge gauge coupling. Adding the contributions 
from the strongly interacting KK quark sector has therefore a much larger impact than 
the associated increase in the effective number of degrees of freedom g e ff- As a result, a e jf 
increases and Qh 2 decreases, as shown in Figs. EH andEb. In contrast, in the case of Z\ LKP, 
the self-annihilation cross sections by themselves are already larger, due to the larger value 
of the weak gauge coupling. The gain from the addition of the KK quark coannihilation 
processes is more than compensated by the associated increase in the effective number of 
degrees of freedom g e ff. As a result, in this case a e // decreases and flh 2 increases, as shown 
in Figs. [3b, and [3b. 

In conclusion, we should mention that the KK Higgs boson Hi in principle can also be 
a potential dark matter candidate. The calculation of its relic density is somewhat more 
model-dependent and we do not consider it here. 



C. Elastic Scattering Cross Sections 



The elastic scattering of the LKP on a nucleon is described by the diagrams depicted 



in Fig. HI For 71 LKP, the corresponding results can be found in 26|, [34| . We follow the 
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FIG. 4: Xree-Ievei diagrams for the elastic scattering of 71 LKP with quarks. The diagrams for 
the case of Z\ LKP are similar. 



computation done in 26j| 7 . The spin-independent cross section is given by 



^"scalar 



Air (m 7l + m T ) : 



[Zf p + (A - Z)f n ] 



(15) 



where my is the mass of the target nucleus, Z and A are respectively the nuclear charge 
and atomic number, while 



u,d,s 



m q 

u,d,s * 



(16) 



and similarly for / n . In eq. ( fl6|) m p (m n ) stands for the proton (neutron) mass. For the 
nucleon matrix elements we take = 0.020±0.004, /£ = 0.026±0.005, = 0.014±0.003, 
/« d = 0.036 ± 0.008, and = 0.118 ± 0.062 Q. The numerical coefficients (3 q and 7 9 in 
eq. ffTBT) are defined as 8 



COS 2 6> w 

„2 



(m 2 i 



m 2 ) 2 



m 2 — m 2 , 

T 1 9l 



COS 2 



K 2 



+ m 2 ! 
(m 2 -m 2 ) 



+ {L^R) 



for a = 0, 



l q = m q 



2 cos 2 ' 



R)\(17) 

(18) 
(19) 



7 The precise calculation of the heavy quark contribution to the processes of Fig. [4] is rather involved - the 
heavy flavors contribute only at the loop level, through the gluon content of the nucleon. In the absence of 
an exact calculation of these effects in the literature, we choose to conservatively ignore the heavy flavor 
contributions altogether, as was done in [26 1. 



Ref. 



34| contains a typo in the overall sign of the coefficient (3 q , which was denoted there as S q . 
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where e is the electric charge, 9w is the Weinberg angle, m q i (m q i) is the mass of an 
S , C/(2)jy-doublet (5'L r (2)vK-singlet) KK quark, and a is the mixing angle in the KK quark 
mass matrix given by sin 2a = 1m q j {Tn q i +m q -Q. Eq. (I17p includes the mixing effect between 
two KK quarks and eq. (|18|) is obtained in the limit when a = 0. This mixing effect gives 
a minor correction to the cross section (at a few percent level) and we do not include it in 



our figures for 5D. However it is important to keep it in the 6D shown in Ref. [251 ] . 

Our convention for the SM hypercharge is Yi = Qi — where Qi (Izi) is the electric charge 
(weak isospin) of particle i. E q in eq. ([18]) is the energy of a bound quark and is rather 
ill-defined. In evaluating eq. ( 1T61) . we conservatively replace E q by the current 9 mass m q . As 
alluded to earlier, in eq. (1181) we only sum over light quark flavors, thus neglecting couplings 
to gluons mediated by heavy quark loops. Note that the two contributions (|18|) and ffTDl to 
the scalar interactions interfere constructively: even with extremely heavy KK quark masses 
(large A qi ), there is an inescapable lower bound on the scalar cross section for a given Higgs 
mass, since the Higgs contribution from eq. ffT9l) scales with the SM Higgs mass and not 
the KK quark masses. 

The analogous results for the case of Z\ LKP can now be obtained from the above formulas 
by simple replacements: m 7l — > m Zl , Y qL — > | and Y qR — > 0, since Zi is mostly the neutral 
SU{2) W gauge boson, which has no interactions with the S , [/(2) M /-singlet KK quarks (or 
equivalently, the right-handed SM quarks). In addition, one should replace co ^ 9w — > sin g w 
to account for the different gauge coupling constant. 

Theoretical predictions for the spin-independent LKP-nucleon elastic scattering cross 
sections are shown in Fig. [5] for different fixed values of the KK quark - LKP mass splitting 
A gi , and for two different LKPs: (a) 71 and (b) Z\. In both cases the cross sections decrease 
as a function of LKP mass. This is due to the inverse scaling of the KK quark exchange 
contributions (fISl) with the KK mass scale. Comparing Fig. [5^ to [5b, we notice that the 
scalar cross section for Z\ is more than one order of magnitude larger than the scalar cross 
section for 71 of the same mass. This is mostly due to the larger SU(2)w gauge coupling. 
Notice that even when the KK quarks are very heavy, there is still a reasonable cross section, 
which is due to the Higgs mediated contribution ffT9l) . Perhaps the most noteworthy feature 



The actual choice of the value for m q is inconsequential since the m q factor in eqs. (fT5|) and (fTT)| cancels 
against the m q factor in the denominator of eq. (|16[) . 
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FIG. 5: Spin-independent elastic scattering cross sections as a function of LKP mass for (a) 71 and 
(b) Z\. The individual curves are labelled by the value of the mass splitting A qi . The SM Higgs 
mass rrih is fixed to 120 GeV. 

of Figs. [5^ and [5b is the significant enhancement of the direct detection signals at small A 9l , 
often by several orders of magnitude. This greatly enhances the prospects for detecting KK 
dark matter, if the mass spectrum turns out to be rather degenerate. 
The spin- dependent cross section is given by 



1 



rrin 



a. 



spm 



6ir (m 71 + %) 



;Jn{Jn + 1, 



-1 2 



u.d.t 



(20) 



where a q and X q are 



2e 2 



a 



cos 2 6 



2 

m \ 



+ {L^R) 



71 



A,= AP(S p )/J N + A n q (S n )/J N 



(21) 
(22) 



Here Jn is the nuclear spin operator. A^ ,ri is given by (p, n\S^\p, n) = AP' n Sg n and is 
the fraction of the nucleon spin carried by the quark q. We use = A^ = 0.78 ± 0.02, 
A^ = A™ = -0.48 ± 0.02, and Af = A™ = -0.15 ± 0.02 [62j. (S p<n )/J N = (N\S p>n \N) / J N is 
the fraction of the total nuclear spin Jn that is carried by the spin of protons or neutrons. 
For scattering off protons and neutrons, \ q reduces to A p q and A^, respectively. 
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Following 63j, we can rewrite eq. (I20I) in the form 



spin 



-G|/i 2 ^-^ (a p (S p > + a n <^» 2 , (23) 

71 J N 



where Gf is the Fermi constant and 



= m T m 7l 
my + m 7l 



is the reduced mass, while the coefficients a p and a n are given by 



8v^3G^m 7l 

^ V 

4V3G F cos 2 w ^ s 



F 2 



m 2 ^ — m 7i 



Ar • (25) 



The main advantage of introducing the parameters a p and a n is that they encode all the 
theoretical model-dependence, thus allowing different experiments to compare their sensi- 
tivities in a rather model-independent way. From eqs. (12311241) it is clear that for any given 
target, the spin- dependent scattering rate depends on only three parameters: m 71 , a p and 
a n . Notice that in our setup there are only two relevant model parameters: m^xp and A qi , 
therefore we will have a certain correlation between a p and a n , depending on the nature of 
the LKP 10 . 

In Fig. [6] we show our result for the spin-dependent LKP elastic scattering cross sections 
off protons and neutrons for the case of (a) 71 and (b) Zi, for different mass splittings A qi . 
The red solid curves are the LKP-proton cross sections and the blue dotted curves are the 
LKP-neutron cross sections. All curves exhibit the same general trends as the corresponding 
spin-independent results from Fig. [5j the cross sections decrease with the KK mass scale, 
and are enhanced for small mass splittings A qi . One peculiar feature is that the proton and 



10 



In introducing the parameters a p and a n we have followed the convention of Ref. [631 ]. We should 



34], where dp n WHS defined aS dp^ri 



alert the reader that a different convention was used in Ref. 
Y^, u 4 s \Yq L + Yq R ] Ag' n i so that a p and a n are pure numerical factors, e.g. a n = — 0.f39f67 and 
a p = 0.280833 for ^1 LKP, and a n = a p = 0.0375 for Z\ LKP. However, this factorization can only 
be done in the special case of Aq x = A qi , and furthermore, the theoretical model dependence (through 
the KK quark masses) creeps back explicitly in the expression for the spin-dependent cross section. 
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FIG. 6: The spin-dependent elastic scattering cross sections as a function of LKP mass for (a) 
71 and (b) Z\ at various A 9l . The red solid curves are the LKP-proton cross sections and the 
blue dotted curves are the LKP-neutron cross sections. The LKP-proton and LKP-neutron cross 
sections are identical for Z\ . For 71 the proton cross section is approximately 4 times larger than 
the neutron cross section, for the same values of the LKP mass m 71 and mass splitting A 9l . 

neutron spin-dependent cross sections are equal in the case of Z±, as seen in Fig. [6b. This 
is an exact statement, which is due to the fact that Zi does not particularly discriminate 
between the different quark flavors in the nucleon - it couples with equal strength to both 
up- and down-type (left-handed) quarks. On the other hand, 71 couples differently to u and 
d, because of the different hypercharges of the right-handed quarks. As a result, the cross 
sections on protons and neutrons differ in the case of 71, as seen in Fig. [6^. Interestingly, for 
a given LKP mass m 7l and mass splitting A qi , the proton cross section in Fig. [6^ is larger 
than the neutron cross section by about a factor of 4, which is due to a numerical coincidence 
involving the values of the quark hypercharges and the A^ ,n parameters. 11 Because of this 
simple scaling, for a given LKP mass m 7l , the proton cross section at a certain A 5l coincides 
with the neutron cross section for half the mass splitting (A qi /2) since to leading order both 
the proton and the neutron cross sections are proportional to (A 9l )~ 2 . 

We shall now review the corresponding results for the case of two universal extra dimen- 
sions. The sum (3 q + 7 g for the spinless photon {^ H ) LKP was computed in [25| (note that 

11 This can be simply understood in terms of the relative scaling of the a p and a n parameters introduced in 
Eq. (|2U)) . In the case of 71, they differ by a factor of -2, while for the case of Z\ they are the same. 
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FIG. 7: The spin-independent elastic scattering cross sections as a function of LKP mass for (a) 
jh in 6D UED and (b) Hi in 5D UED. In the case of 'Jh LKP (panel "a"), we present results for 
several values of A qi as shown. In the case of Hi (panel "b"), we only show the contribution from 
the SM Higgs exchange, assuming heavy KK quarks. 



here we are using a different convention for the hypercharges Y{) 
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Pq +lq 



cos 2 9 W 
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quJ 
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m? - (m q + m jH ) 2 m? - (m 



+ 



m q 
2<j 



,(26) 



where m lH is the mass of the spinless photon, m qi is the (common) mass of the SU(2)w- 
doublet and SU (2) ^-singlet KK quarks, while m q is the corresponding SM quark mass. 

Using Eqn. (fl5|) . we obtain the spin- independent elastic scattering cross section for 7^ 
as shown in Fig. [7^. The different curves are labelled by the assumed fixed value of A qi , 
and are plotted versus the LKP mass m lH . We see that the size of the jh signal is about 
the same order as the 71 cross sections from Fig. [5h- On the other hand, the relic density 
constraint would single out somewhat different regions for m 7l and m lH . The annihilation 
cross section for ju is smaller than that of 71 [25j, and correspondingly, lower 7^ masses 
would be preferred, with enhanced prospects for direct detection 12 . Notice that there is no 



12 One should keep in mind that mass splittings A gi as small as those shown in the plot mandate the 
inclusion of coannihilation processes with KK quarks in 6D UED, as we did previously in the case of 71 
and Z\ in 5D UED (see Fig. [1} . Unfortunately the computation of coannihilations in 6D UED does not 
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spin-dependent cross section for 7# since it is a scalar particle. 

In conclusion of this section, we shall briefly discuss the scenario of KK Higgs (Hi) 
LKP. Just like 7#, Hi is a scalar and does not have spin- dependent interactions. Its spin- 
independent elastic scattering cross section can be readily computed following the procedure 
outlined earlier in this section and in the appendix. In this case, the KK quark exchange 
diagrams are also Yukawa suppressed, and the dominant among them is the s KK quark 
contribution. As in the 71 LKP case, the diagrams with KK quark exchange and SM 
Higgs exchange interfere constructively. Therefore, the SM Higgs exchange diagram by itself 
provides a conservative lower bound on the elastic scattering cross section, independent of 
the other details of the KK spectrum, and in particular, the KK quark masses. This absolute 
minimum of the cross section is plotted in Fig. [7b as a function of the LKP mass mn 1 ■ It 
is worth mentioning that this result is completely independent of the SM Higgs mass mn- 
The contribution corresponding to (|T9|) is given by 

3 e 2 m q 

7<? = Lin 2 a ' ( 2 ') 

4 Sill U\y "lyy 



where mjy is the mass of the W boson. The coupling of the KK Higgs to the SM Higgs 

A- 



boson is the same as the triple Higgs coupling of the SM, which is proportional to m 2 . This 



m\ dependence is exactly cancelled by the m^ 2 dependence of the SM Higgs propagator 
in the non-relativistic limit (see Eqn. ( TL9i) ). Therefore the final cross section is indeed 
independent of the SM Higgs mass, and this fact remains true regardless of the values of the 
KK quark masses. 

III. DIRECT WIMP DETECTION AND EXPERIMENTS 

The detailed distribution of dark matter in our galaxy, and in particular in the local 
neighborhood, is not well constrained by current observations and high-resolution simu- 
lations. The standard assumption for its distribution is a cored, non-rotating isothermal 
spherical halo with a Maxwell-Boltzmann velocity distribution with a mean of 220km/s, 



and escape velocity from the galactic halo of 544km/s 64j. For the local density of dark 



exist in the literature. In analogy with the case of 71 LKP in 5D UED, we expect the coannihilation 
effects to increase the preferred range of m 1H . 



24 



matter particles we assume p=0.3GeV/cm 3 65]. 



661 ] con- 



The WIMP interaction signature in ultra-low-background terrestrial detectors 
sists of nuclear recoils. Direct detection experiments attempt to measure the small 
(<100keV) energy deposited when a WIMP scatters from a nucleus in the target medium. 
The recoil energy of the scattered nucleus is transformed into a measurable signal, such 
as charge, scintillation light or lattice excitations, and at least one of the above quanti- 
ties can be detected. Observing two signals simultaneously yields a powerful discrimina- 
tion against background events, which are mostly interactions with electrons as opposed 
to WIMPs and neutrons, which scatter from nuclei. The WIMP interaction takes place in 
the non-relativistic limit, therefore the total cross section can be expressed as the sum of 
a spin-independent (SI) part (see Eqn. ffT5l) ). a coherent scattering with the whole nucleus, 
and of a spin-dependent (SD) part (see Eqn. (1201) ). which describes the coupling to the total 
nuclear spin [671 ] . 

Neutrons with energies in the MeV range can elastically scatter from nuclei and mimic 
a WIMP signal. Two methods are used to discriminate against the residual neutron back- 
ground, which comes from (a,n)- and fission-reactions in materials and from interactions of 
cosmic muons with the rock and experimental shields. First, the SI WIMP-nucleus cross 
section is proportional to the atomic mass-squared of the nucleus, making the expected 
total WIMP interaction rate material dependent. Second, the mean free paths of WIMPs 
and MeV neutrons are exceedingly different (10 10 m versus 8 cm in a typical WIMP target), 
allowing to directly constrain the neutron background from the ratio of observed single to 
multiple interaction events. 

The experimental upper bounds of the SI cross section from direct detection experiments 
are WIMP-type independent and thus will not change if we consider different WIMP can- 
didates. Similarly, the SD cross section limits can also be reinterpreted for various DM 
candidates. The only exception is a spin zero WIMP, such as 7# in 6D UED, which does 
not have an axial-vector coupling with nuclei, hence no SD interaction is expected. We will 
extensively discuss the model dependence of the SD cross section in the next section. 

In this study, we choose four direct detection experiments which demonstrated best ex- 
perimental sensitivity to-date in various parts of the WIMP search parameter space. The 



CDMS experiment sets the best SI upper bound above a WIMP mass of 42 GeV 



68(| . while 



XENON10 gives the most stringent upper bound on WIMP-neutron SD couplings 7_0( and 
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Total 


F/nprP'v ranP"P TjOPation F»pf 


CDMS II 


Ge(Si) 


4.75 kg (1.1kg 


) lOkeV- lOOkeV Soudan, USA [68] 


XENON10 


Xe 


15 kg 


4.5 keV - 27keV Gran Sasso, Italy [69] 


KIMS 


Csl 


34.8 kg 


3keV - llkeV Yangyang, Korea [71] 


COUPP 


CF 3 I 


1.5 kg 


5 keV - Fermilab, USA [72] 



TABLE I: Direct WIMP detection experiments considered in this study. 
SI couplings below 42 GeV 



69|. The KIMS |7l| and COUPP [72] experiments show the 
best sensitivity for SD WIMP-proton couplings. As we shall see in the following section, the 
combined study of all four experiments strongly constrains the SD proton-neutron mixed 
coupling parameter space (the so called a p -a n parameter space, where a p and a n are the dark 
matter particle's couplings to protons and neutrons, respectively, see eq. (1231)). 

Table [I] summarizes the relevant characteristics of the four experiments such as target 
material, total mass, energy range considered for the WIMP search, and location. In this 
paper we either calculated the LKP limits based on published data (XENON10), or we 
obtained the data points for the cross section upper bounds from the collaboration (CDMS, 
KIMS and COUPP). 

The CDMS experiment [68] is operated in the Soudan Underground Laboratory, USA. 
It uses advanced Z(depth)-sensitive Ionization and Phonon (ZIP) detectors, which simul- 
taneously measure the ionization and athermal phonon signals after a particle interacts in 
the crystal. The ZIP detectors provide excellent event-by-event discrimination of nuclear 
recoils from the dominant background of electron recoils. The most stringent limits on spin- 
independent couplings with nucleons above a WIMP mass of 42 GeV comes from the first two 
CDMS-II five tower runs with a raw exposure of 397.8 kg-days in germanium. The null ob- 
servation of a WIMP signal sets a WIMP-nucleon cross section upper bound of 6.6xl0 _8 pb 
(for a 60 GeV WIMP mass) and of 4.6 x 10 -8 pb when the results are combined with previous 
CDMS results. 

The SuperCDMS project [l^Ql is a three-phase proposal to utilize CDMS-style detectors 
with target masses growing from 25 kg to 150 kg and up to 1 ton, with the aim of reaching 
a final sensitivity of 3xl0 _11 pb by mid 2015. This goal will be realized by developing 
improved detectors and analysis techniques, and concomitantly reducing the intrinsic surface 
contamination of the crystals. 



The XENON10 collaboration 



69J operated a 15 kg active mass, dual-phase (liquid and 
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gas) xenon time projection chamber in the Gran Sasso Underground Laboratory (LNGS), 
in WIMP search mode from August 2006 to February 2007. XENON10 uses two arrays of 
UV-sensitive photomultipliers (PMTs) to detect the prompt and proportional light signals 
induced by particles interacting in the sensitive liquid xenon (LXe) volume. The 3D position 
sensitivity, the self-shielding of LXe and the prompt versus proportional light ratio are 
the most important background rejection features. The first results, using ~136kg-days 
exposure after cuts, demonstrated that LXe can be used for stable, homogeneous, large scale 
dark matter detectors, providing excellent position resolution and discrimination against the 
electron recoil background. The derived upper bound on SI cross sections on nucleons is 
4.5 xl0~ 8 pb for a WIMP mass of 30 GeV. Since natural Xe contains 129 Xe (26.4%) and 
131 Xe (21.2%) isotopes, each of these having an unpaired neutron, the XENON10 results 
substantially constrain the SD WIMP-nucleon cross section. We calculated the XENON10 
SD LKP-neutron and LKP-proton upper bounds based on the observation of 10 events, 
without any background subtraction 70]. The next phase, XENON100, will operate a total 
of 170 kg (70 kg fiducial) of xenon, viewed by 242 PMTs, in a dual-phase TPC in an improved 
XENON10 shield at the Gran Sasso Laboratory. While the fiducial mass is increased by 
more than a factor of 10, the background will be lower by about a factor of 100 (through 
careful selection of ultra-low background materials, the placing of cryogenic devices and 
high- voltage feed-throughs outside of the shield and by using 100 kg of active LXe shield) 
compared to XENON10. XENON100 is currently being commissioned at LNGS, the aim is 
to start the first science run in fall 2008, probing WIMP-nucleon SI cross sections down to 
~10~ 9 pb. 

The Korea Invisible Mass Search (KIMS) experiment 7lJ is located at the Yangyang Un- 
derground Laboratory, Korea. The collaboration has operated four low-background CsI(Tl) 
crystals, each viewed by two photomultipliers, for a total exposure of 3409kg-days. Both 
133 Cs and 127 I are sensitive to the spin-dependent interaction of WIMPs with nuclei. KIMS 
detects the scintillation light after a particle interacts in one of the crystals, kept stably at 
(0±0.1)°C. The pulse shape discrimination technique, using the time distribution of the 
signal, allows to statistically separate nuclear recoils from the electron recoil background. 
The KIMS results are consistent with a null observation of a WIMP signal, yielding the 
best limits on SD WIMP-proton couplings for a WIMP mass above 30 GeV. Specifically, the 
upper bound for a WIMP mass of 80 GeV is 1.7xl0 -1 pb. 
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FIG. 8: Current and projected experimental limits on the SI LKP-nucleon-scattering cross section 
together with the theoretically expected 71 (blue shaded), Z\ (yellow shaded) and jh (green- 
shaded) LKP regions. The boundaries of the LKP regions are selected for 0.01 < A < 0.5 while 
the Higgs mass is hxed to 120 GeV. The solid lines are the current experimental upper bounds 
(90% C.L) from the CDMS (blue) and XENON10 (red) experiments. The dashed lines are expected 
sensitivities for the SuperCDMS 25 kg (blue) and XENON100 (red) experiments, which will be 
operated in the near future. The dotted line is the expected sensitivity for a ton-scale detector. 

The Chicagoland Observatory for Underground Particle Physics (COUPP) experi- 
ment [72I is operated at Fermilab, USA. The experiment has revived the bubble chamber 
technique for direct WIMP searches. The superheated liquid can be tuned such that the 
detector responds only to keV nuclear recoils, being fully insensitive to minimum ionizing 
particles. A 1.5 kg chamber of superheated CF 3 I has been operated for a total exposure of 
250kg-days. The presence of fluorine and iodine in the target makes COUPP sensitive to 
both SD and SI WIMP-nucleon couplings. The production of bubbles is monitored optically 
and via sound emission, reaching a reconstructed 3D spatial resolution of ~1 mm. It allows 
to reject boundary-events and to identify multiple neutron interactions. The most recent 
COUPP results set the most sensitive limit on SD WIMP-proton cross sections for a WIMP 
mass below 30 GeV. As an example, the upper bound on the SD coupling is 2.7xl0 _1 pb at 
a WIMP mass of 40 GeV. 

In Fig.[8]we show the current CDMS and XENON10 upper bounds for the SI cross section 
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FIG. 9: Current experimental limits on the SD (a) neutron- and (b) proton-scattering cross section 
together with the predicted SD LKP-neutron (proton) cross sections for 71 (blue shaded) and Z\ 
(yellow shaded). The solid curves for each plot are the upper bounds (90% C.L.) from the COUPP 
(green), KIMS (black) and XENON10 (red) experiments. The dotted line shows the expected 
sensitivity for a ton-scale detector which is obtained from Fig. [8] by a proper rescaling considering 
a Xenon detector. The 71 and Z\ LKP areas are drawn with the same A qi convention as shown in 

Fig. m 



together with projected sensitivities for SuperCDMS 25 kg, XENON100 and for a ton-scale 
detector. The LKP boundaries for 71, Z\ and 7^ as dark matter candidates are also shown, 
for a wide range of mass splittings (0.01 < A 91 < 0.5) and a fixed Higgs mass nth of 120 GeV. 
The small mass splitting regions are excluded up to a mass of about 600 GeV, 900 GeV and 
700 GeV for 7 1; Zi and 7#, respectively. For large mass splittings of A qi = 0.5, only masses 
below about 100 GeV can be probed. Future ton-scale direct detection experiments should 
cover most of the interesting LKP parameter space. 

In Fig. M, we show the SD cross section limits for both (a) pure neutron and (b) pure 
proton couplings for three experiments together with the theoretical predictions for 71 and 
Z\ for a range of mass splittings (0.01 < A qi < 0.5). The most stringent SD pure neu- 
tron upper bound is set by the XENON10 experiment, while the best SD cross section for 
pure proton couplings in the region of interesting LKP masses (> 500 GeV) comes from the 
KIMS experiment. As explained in the previous section, the theoretical 71 and Z\ regions 
are overlapping for pure neutron couplings, while for pure proton coupling these can be 
distinguished for a given mass splitting A gi . 

In the following section we investigate the details of the LKP specific parameter spaces. 
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IV. LIMITS ON KALUZA-KLEIN DARK MATTER 



In the previous sections we introduced the different dark matter candidates in UED 
models: KK gauge bosons (71 and Zi) and KK scalars (7^ and Hi). On the theoretical side, 
we discussed the calculation of their relic densities and elastic scattering cross sections. On 
the experimental side, we described the different types of experiments which are sensitive to 
KK dark matter. We shall now combine our theoretical predictions with the current /future 
measurements discussed earlier. Where applicable, we shall also include constraints from 
high energy collider experiments. We shall be particularly interested in the region of small 
mass splittings A qi , which is problematic for collider searches, but promising for direct 
detection. We will concentrate on KK gauge boson dark matter (both 71 and Zi), whose relic 



23, 



13 



density can be reliably calculated, including all relevant coannihilation processes 

In Fig. [10] we present a combination of results for the case of (a) 71 and (b) Zi LKP in 
5D UED. As we emphasized earlier, the two most relevant parameters are the LKP mass 
(m 7l or mz 1 , correspondingly) and the mass splitting A 9l between the LKP and the KK 
quarks. We therefore take both of these parameters as free and do not assume the MUED 
relation among them. For simplicity, we assume that the SU (2)v^-doublet KK quarks and the 
SU (2) w-sin.g\et KK quarks are degenerate, so that there is a single mass splitting parameter 
which we have been calling A qi . However, this assumption is only made for convenience, and 
does not represent a fundamental limitation - all of our results can be readily generalized 
for different KK quark mass splittings (i.e. several individual A parameters). The masses of 
the remaining KK particles in the spectrum are fixed as in Fig. [TJ in the case of 71 LKP, we 
use the MUED spectrum, while in the case of Z\ LKP, we take the gluon and the remaining 
particles to be respectively 20% and 10% heavier than the Zi. This choice is only made for 
definiteness, and does not carry a big impact on the validity of our results, as long as the 
remaining particles are sufficiently heavy so that they do not participate in coannihilation 
processes. 

In the so defined parameter plane, in Fig. [TU] we superimpose the limit on the spin- 
independent elastic scattering cross section, the limit on the relic abundance and the LHC 



reach in the four leptons plus missing energy (AC+Efx) channel which has been studied in [45 ]. 



13 The scalar candidates (jh and Hi) do not have spin-dependent interactions anyway, while their spin- 
independent scattering rates are similar to the examples we consider, as discussed in Sec. Ill CI 
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FIG. 10: Combined plot of the direct detection limit on the spin-independent cross section, the 
limit from the relic abundance and the LHC reach for (a) 71 and (b) Z\, in the parameter plane 
of the LKP mass and the mass splitting A qi . The remaining KK masses have been fixed as in 
Fig. [I] and the SM Higgs mass is = 120 GeV. The black solid line accounts for all of the 
dark matter (100%) and the two black dotted lines show 10% and 1%, respectively. The green 
band shows the WMAP range, 0.1037 < ttcDHih 2 < 0.1161. The blue (red) solid line labelled 
by CDMS (XENON 10) shows the current limit of the experiment whereas the dashed and dotted 
lines represent projected limits of future experiments as shown in Fig. In the case of 71 LKP, 
a ton-scale experiment will rule out most of the parameter space while there is little parameter 
space left in the case of Z\ LKP. The yellow region in the case of 71 LKP shows parameter space 
that could be covered by the collider search in the Ai + Fj>T channel at the LHC with a luminosity 
of 100 fb~ l 0]. 



This signature results from the pair production (direct or indirect) of S'[/(2)iy-doublet KK 
quarks, which subsequently decay to Z\S and jets. The leptons (electrons or muons) arise 
from the Z\ — > £ + £~ , ~fi decay, whose branching fraction is approximately 1/3 [45] . Requiring 



a 5<t excess at a luminosity of 100 fb x , the LHC reach extends up to R 1 



m 7l ~ 



1.5 TeV, 



which is shown as the right-most boundary of the (yellow) shaded region in Fig. [TOh . The 
slope of that boundary is due to the fact that as A qi increases, so do the KK quark masses, 
and their production cross sections are correspondingly getting suppressed, diminishing the 
reach. We account for the loss in cross section according to the results from Ref. 75], 
assuming also that, as expected, the level-2 KK particles are about two times heavier than 
those at level 1. Points which are well inside the (yellow) shaded region, of course, would be 
discovered much earlier at the LHC. Notice, however, that the LHC reach in this channel 
completely disappears for A gi less than about 8%. This is where the KK quarks become 
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lighter than the Z\ (recall that in Fig. [TOa mz 1 was fixed according to the MUED spectrum) 
and the q\ — > Z\ decays are turned off. Instead, the KK quarks all decay directly to the 71 
LKP and (relatively soft) jets, presenting a monumental challenge for an LHC discovery. So 
far there have been no studies of the collider phenomenology of a Z\ LKP scenario, but it 
appears to be extremely challenging, especially if the KK quarks are light and decay directly 
to the LKP. This is why there is no LHC reach shown in Fig. [TOb . In conclusion of our 
discussion of the collider reaches exhibited in Fig. [101 we draw attention once again to the 
lack of sensitivity at small A q , : such small mass splittings are quite problematic for collider 



searches (see, for example, [76|, |77|] for an analogous situation in supersymmetry) . 

In Fig. [TU] we contrast the LHC reach with the relic density constraints and with the 
sensitivity of direct detection experiments. To this end we convert our results from Figs. [JJ 



and [8] into the m^KP-A qx plane shown in Fig. [TUl The green s 
100% represents 2a WMAP band, 0.1037 < Vt CDM h 2 < 0.1161 



laded region labelled by 
and the black solid line 



inside this band is the central value flcDMh 2 = 0.1099. The region above and to the right 
of this band is ruled out since UED would then predict too much dark matter. The green- 
shaded region is where KK dark matter is sufficient to explain all of the dark matter in the 
universe, while in the remaining region to the left of the green band the LKP can make up 
only a fraction of the dark matter in the universe. We have indicated with the black dotted 
contours the parameter region where the LKP would contribute only 10% and 1% to the 
total dark matter budget. Finally, the solid (CDMS in blue and XENON10 in red) lines 
show the current direct detection limits, while the dotted and dashed lines show projected 
sensitivities for future experiments (for details, refer back to Sec. IIIII) 14 . 

Fig. [10] demonstrates the complementarity between the three different types of probes 
which we are considering. First, the parameter space region at very large miKP is inconsis- 
tent with cosmology - if the dark matter WIMP is too heavy, its relic density is too large. 
The exact numerical bound on the LKP mass may vary, depending on the particle nature 
of the WIMP (compare Fig. [TUh to Fig. [TUb) and the presence or absence of coannihilations 
(compare the uilkp bound at small A qi to the bound at large A qi ). Nevertheless, we can 



14 Here and in the rest of the paper, when presenting experimental limits in an underdense or an overdense 
parameter space region, we do not rescale the expected direct detection rates with the calculated relic 
density. The latter is much more model-dependent, e.g. the mismatch with the WMAP value may be 
fixed by non-standard cosmological evolution, having no effect on the rest of our analysis. 
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FIG. 11: Limit on the SM Higgs mass and the LKP mass for (a) 71 and (b) Z\ for a given 
A qi = 0.1. The horizontal black solid line is the lower bound on mh from LEP-II (90% C.L. at top 
quark mass 173 GeV). The diagonal black solid line delineates the region disfavored by precision 
data. The solid curves are the current (90% C.L) limits from CDMS (in blue) and XENON10 (in 
red). The dashed curves (SuperCDMS 25kg and XENON100) and dotted line (ton-scale detector) 
are the projected sensitivities for the future experiments. 



see that, in general, cosmology does provide an upper limit on the WIMP mass. On the 
other hand, colliders are sensitive to the region of relatively large mass splittings A 9l , while 
direct detection experiments are at their best at small A qi and small vn^p. The relevant 
parameter space is therefore getting squeezed from opposite directions and is bound to be 
covered eventually. This is already seen in the case of 71 LKP from Fig. [TUk : the future ex- 
periments push up the current limit almost to the WMAP band. Unfortunately in the case 
of Z\ LKP the available parameter space is larger and will not be closed with the currently 
envisioned experiments alone. However, one should keep in mind that detailed LHC studies 
for that scenario are still lacking. 

While previously we already argued that mix p and A 91 are the most relevant parameters 
for UED dark matter phenomenology, for completeness we also investigate the dependence on 
the SM Higgs mass rrih, which is currently still unknown. In Fig.[Tl]we therefore translate the 
information from Fig.[8]into the mLKP-mh plane, for a given fixed KK mass splitting A = 0.1 
now taking the Higgs mass as a free parameter. In each panel, the horizontal black solid 
lines mark the current Higgs mass bound of 114 GeV while the diagonal black solid lines show 
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the indirect limit from the oblique corrections in this model 59|. 15 For low m/,, the limit on 
the LKP mass (or equivalently, the compactification scale) is ttilkp ~ R~ l > 600 GeV (for 
m t = 173 GeV), but it gets weaker for larger m/,, so that ttllkp values as low as 300 GeV are 
still allowed if the SM Higgs boson is very heavy [60J. In Fig. [TT]we also show the current 
(solid lines) limits from CDMS (in blue) and XENON10 (in red), their projected near-future 
sensitivities, SuperCDMS 25 kg and XENON100 (dashed lines), and the projected sensitivity 
of a ton-scale detector (dotted line). The shape of these contours is easy to understand. At 
large m^, the Higgs exchange diagram in Fig. H] decouples, the elastic scattering rate becomes 
independent of nth and the direct detection experimental sensitivity is only a function of 
ttilkp (since A qi is held fixed). In the other extreme, at small rrih, the Higgs exchange 
diagram dominates, and the sensitivity now depends on both rrih and ttilkp- Unfortunately, 
for A = 0.1 the current direct detection bounds do not extend into the interesting parameter 
space region, but future experiments will eventually start probing the large rrih corner of the 
allowed parameter space. On the positive side, one important lesson from Fig. [TT] is that 
the rrih dependence starts showing up only at very low values of rrih, which have already 
been ruled out by the Higgs searches at colliders. This observation confirms that when it 
comes to interpreting existing and future experimental limits on WIMPs in terms of model 
parameters, ttllkp and A 9l are indeed the primary parameters, while rrih plays a rather 
secondary role. 

We remind the reader that the LHC will be able to probe all of the parameter space 
shown in Fig. [TTh through the 4£ + E^t signature, while the discovery of UED in Fig. [TTb 
appears quite problematic. Of course, the SM Higgs boson will be discovered in both cases, 
for the full range of rrih masses shown. 

We now turn to a discussion of the corresponding spin-dependent elastic scattering cross 
sections, which also exhibit an enhancement at small A qi , as shown in Fig. [6j Similar to 
Fig. [lOl in Fig. [12] we combine existing limits from three different experiments (XENON10, 
KIMS and COUPP) in the rriLKP-A qi plane. Panel (a) (panel (b)) shows the constraints 
from the WIMP-neutron (WIMP-proton) SD cross sections. The rest of the KK spectrum 



15 One should keep in mind that the latter have been calculated only for the case of LKP, and only 
within the framework of minimal UED. The line shown in Fig. II lb is therefore only for illustration. 
Furthermore, the 71 calculation itself may be subject to modifications in the more general scenarios which 
we are considering here. 
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FIG. 12: Experimental upper bounds (90% C.L.) on the spin-dependent elastic scattering cross 
sections on (a) neutrons and (b) protons in the niLKP-A qi plane. The solid (dashed) curves are 
limits on 71 (Z\) for each experiment. Shaded regions and dotted lines are defined in the same way 
as in Fig. [7D1 The depicted LHC reach (yellow shaded region) applies only to the case of 71 LKP. 



has been fixed as in Figs. [I] and [OH and m h = 120 GeV. The solid (dashed) curves are limits 
on 7x (Zi) for each experiment. The constraints from LHC and WMAP on the m^p-A,^ 
parameter space are the same as in Fig. [TUJ 

By comparing Figs. [TD] and [T5] we see that, as expected, the parameter space constraints 
from SI interactions are stronger than those from SD interactions. For example, in perhaps 
the most interesting range of LKP masses from 300 GeV to 1 TeV, the SI limits on A qi in 
Fig. [TU1 range from a few times 10~ 2 down to a few times 10~ 3 . On the other hand, the SD 
bounds on A gi for the same range of mixp are about an order of magnitude smaller (i.e. 
weaker). We also notice that the constraints for 71 LKP are stronger than for Z\ LKP. This 
can be easily understood by comparing Fig. and Fig. Eb: for the same LKP mass and 
KK mass splitting, the 71 SD cross sections are typically larger. 

Fig. [12] also reveals that the experiments rank differently with respect to their SD limits 
on protons and neutrons. For example, KIMS and COUPP are more sensitive to the proton 
cross section, while XENON10 is more sensitive to the neutron cross section. As a result, 
the current best SD limit on protons comes from KIMS, but the current best SD limit on 
neutrons comes from XENON10. Combining all experimental results can give a very good 
constraint on the a p -a n parameter space. Fig. [T5k (Fig. fT3"b) shows combined results for 
TfiLKP = 50 GeV {rriLKP = 500 GeV) in the (model-independent) a p -a n parameter space. 
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FIG. 13: Experimental limits on the a p -a n parameter space for (a) miKP = 50 GeV and (b) 
m LKP = 500 GeV. The contours show limits from XENON10 (red solid line), KIMS (black dotted 
line) and COUPP (green dashed line). The blue near-horizontal bands show the evidence regions 
allowed by DAMA \7a\ . while the green region shows the parameter space allowed by all current 
experimental results. The two straight lines originating from a n = a p = are the theoretical pre- 
dictions for dp and a n in the case of 71 or Z\ LKP in 5D UED. These theory lines are parametrized 
by the value of A qi as indicated by a few representative points. 



The contours show limits from XENON10 (red solid line), KIMS (black dotted line) and 
COUPP (green dashed line). The blue near- horizontal bands show the evidence regions 
allowed by DAMA [78], while the green region shows the parameter space allowed by all 
current experiments. Note that these limits were computed in two different ways. The 
results from KIMS and COUPP are based on the method proposed in 63| whereas those 



781 ] . We believe that the latter is 



from DAMA and XENON10 are calculated as advocated in 
more accurate since limits are computed for all angles in the a p -a n plane separately whereas 
the former solely relies on the limits calculated considering pure coupling to neutrons and 
protons respectively. More details about these calculations can be found in the appendix. 
The two straight lines originating from a n = a p = are the theoretical predictions for a p and 
a n in the case of 71 or Z\ LKP in 5D UED. These theory lines are parametrized by the value 
of A qi as indicated by a few representative points. The feature which is readily apparent 
in Fig. [13] is the orthogonality between the regions allowed by the a p -sensitive experiments 
like KIMS and COUPP, on the one side, and the a n -sensitive experiments like XENON10, 
on the other. This indicates the complementarity of the two groups of experiments: the 
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FIG. 14: Experimental limits (90% C.L.) on the scalar LKP (-y H ) in 6D UED. (a) Lower bound 
of A qi vs m~ fH for rrih =120 GeV. The solid lines are the current experimental lower bounds on 
A qi for a given m lH from CDMS (blue) and XENON10 (red). SuperCDMS 25 kg and XENON100 
projected sensitivities are drawn with dashed lines. The dotted line shows the projected sensitivity 
of a ton-scale experiment, (b) Lower bound of the Higgs mass (m^) as a function of m^ H for a 
fixed A qi =0.1. The WMAP allowed range is the green shaded region. The LEP II lower bound 
on rrih is shown as the black solid line. 

green-shaded region allowed by the combination of all experiments is substantially more 
narrow than the region allowed by each individual experiment. 

In conclusion of this section, we shall also consider KK dark matter candidates in models 
with two universal extra dimensions (6D UED). As mentioned in Sec. Ill Al the novel possi- 
bility here compared to 5D UED is the scalar photon (7^) LKP. As a spin zero particle, it 
has no spin-dependent interactions and can only be detected through its spin-independent 
elastic scattering. Fig. [T4"h (Fig. [T4"b) is the analogue of Fig. [10] (Fig. [IT]) for the case of 7^ 
LKP. In Fig. [T4"h we show lower bounds on A qi versus the mass m 7H of the scalar photon, for 
a fixed Higgs mass (rrih = 120 GeV). The solid lines indicate the current experimental limits 
from CDMS (blue) and XENON10 (red). The dashed lines are the projected sensitivities 
of SuperCDMS 25 kg and XENON100 and the dotted line is the projected sensitivity of 
a ton-scale detector. Since the cosmologically preferred mass range for 7^ is much lower 
(~ 200 GeV before accounting for coannihilations) than for the LKP in 5D UED, the con- 
straints are quite powerful - in particular, the future ton-scale experiments are expected to 
cover most of the interesting mass splitting (A qi ) region. 



37 



In Fig. [Tib we show lower bounds of the Higgs mass as a function of m lH for a fixed 
A qi = 0.1. The WMAP preferred parameter space is marked as the green shaded region, 
while the black solid line is the LEP II lower limit on m^. The contours resemble in shape 
those seen earlier in Fig. [TTJ In particular, we notice that within the LEP II allowed range, 
the Higgs mass does not have a large impact on the direct detection bounds. However, if the 
LHC finds a SM Higgs boson with a mass smaller than ~300 GeV, then the WMAP bound 
would constrain the mass of 7# within a relatively narrow mass ranges at a given mass 
splitting (A gi ). For example in Fig. [Tib , where the fixed mass splitting is A qi = 0.1, the 
corresponding constraint on the mass of 7# would be 180 GeV < m lH < 250 GeV. In fact, 
this conclusion is rather insensitive to the particular choice of A qi . This is due to the fact 
that 7# self-annihilation is helicity-suppressed and gauge boson final states are dominant 
in the WMAP allowed regions. Therefore, Fig. [Tib would look qualitatively similar, if a 
different value of A„, were used. 



V. CONCLUSIONS 



The dark matter puzzle is among the most intriguing questions in particle physics. Its 
origin resides in cosmological observations such as the rotation curves of galaxies, cosmic 
microwave background, gravitational lensing, large scale structure, the mass to luminosity 
ratio and so on. Interestingly, many scenarios of new physics beyond the Standard Model 
provide a stable neutral particle which, in principle, can be produced and observed at col- 
liders. In fact, one of the primary motivations for SUSY has always been the fact that it 
naturally accommodates a WIMP candidate. More recently, we have learned that extra 
dimensional models provide a viable alternative to SUSY dark matter, namely KK dark 
matter. Both of these scenarios have been attracting a lot of attention in terms of collider 
and astrophysical aspects. In this paper we performed a comprehensive phenomenological 
analysis of KK dark matter in universal extra dimensions, extending previous studies by 
considering new LKP candidates [Z\ and Hi). We also revisited the cases of 71 and 7^ 
LKP, focusing on the possibility of a small mass splitting with the KK quarks. All of these 
features can be realized in non-minimal UED scenarios and therefore deserve attention. 

In our analysis we included the relevant theoretical constraints from cosmology (the relic 
density of KK dark matter) and particle physics (low energy precision data). We accounted 
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for all coanniliilation processes in our relic density calculation, focusing on coannihilations 
with KK quarks since they play an important role for direct detection at small mass split- 
tings. 

We then contrasted the sensitivities of the LHC and the different types of direct detection 
experiments, and exhibited their complementarity. We demonstrated that the ttilkp ~ A qi 
parameter space is both convenient and sufficient for a simultaneous discussion of collider and 
direct detection searches. Collider experiments like the LHC and possibly ILC are sensitive 
to the region of relatively low tulkp and sufficiently large A qi . On the other hand, direct 
detection experiments do best at relatively low ttilkp and small A qi . Finally, cosmology 
rules out the region of very large mi,KP- We see that, at least in principle, the combination 
of all three types of constraints has the potential to completely cover the relevant parameter 
space. We showed that with the expected sensitivity of the next generation direct detection 
experiments, the coverage is almost complete in the case of 71 LKP. 

In conclusion, we summarize the main lessons from each of the three main areas in our 
study, and point towards interesting directions for future work. 

• Direct detection of KK dark matter. The direct detection prospects are greatly en- 



hanced when the LKP becomes degenerate with the KK quarks 26]. Therefore, the 
mass splitting A qi is a key parameter which is worth exploiting by the experimen- 
tal collaborations in presenting their limits. The conventional approach is to overlay 
the results of random scans over the full parameter space of a given model over the 
model-independent exclusion curves from Figs. [8] and [91 There are several problems 
with this. First, due to the multitude of parameters being scanned, such scans are 
not sufficiently exhaustive and are almost guaranteed to miss the relevant parts of 
parameter space. For example, suppose we scan a KK quark mass in the range from 
f^LKP up to some maximum value m™ ax . In order to obtain a single parameter space 
point at a given A qi , the number of points sampled along the m qi direction should be 
at least 

_ m™ x - m LK p _ m™ 00 - m LKP 1 1 



m LKP A qi A qi 

In order to probe mass degeneracies of A 9l ~ 0.01, one should therefore generate at 
least 10 2 points along each parameter dimension. It is then clear that random scans 
cannot effectively probe models which have more than a few input parameters. For 
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example, typical SUSY scans often utilize a constrained parameter space which nev- 
ertheless still has on the order of 10-15 parameters. In that case, seeing the effects of 
1% degeneracies would require 10 20 — 10 30 total points - a number which is obviously 
impractical. Furthermore, many of the scanned parameters often have little effect on 
the experimental signals being discussed. We therefore find it much more efficient and 
illuminating to forego such general scans in favor of simple parametrizations which 
would contain only the variables relevant for the experimental search. A simple im- 
plementation of this idea is the m LKP — A ?1 parameter space used in Figs. [T0land[T2l 
Notice that our argument is not limited to UED models - the mass of the dark matter 
particle and the mass splitting are expected to be the two most important parameters 
for dark matter searches in any other theory (such as SUSY, little Higgs, warped extra 
dimensions, etc.). 

• Collider searches for KK dark matter. While beneficial for direct DM detection, mass 
degeneracies are generally problematic for collider searches, since they degrade the 
quality of the discovery signatures. For example, at small A qi the KK quark decay 
products become softer, and the leptonic modes may altogether turn off. Detailed 
LHC studies for these challenging situations are still lacking and are definitely worth 
undertaking. In particular, the LHC phenomenology of Z\ or Hi LKP has never been 
discussed. 

• Theoretical calculations related to KK dark matter. In this paper we have reviewed 
the main ingredients of a complete analysis of KK dark matter phenomenology. We 
emphasized the importance of including the effect of coannihilations in the case of 
small mass splittings which are relevant for experimental searches. However, there are 
still several missing pieces which are needed to improve the accuracy of our predictions. 
While the relic density calculation for 71 and Z\ LKP is on a relatively firm footing, 
since the corresponding coannihilation processes are known, the complete calculation 
of coannihilation effects in case of Hi, 7^ or Zjj LKP is still lacking. Similarly, one 
would like to have available the precise calculation of the heavy flavor contribution 
to the LKP elastic scattering cross section on nucleons. We are hoping that it would 
take less than a WIMP direct detection signal to jumpstart the theoretical efforts in 
this direction. 
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VI. APPENDIX: EVENT RATE CALCULATIONS FOR DIRECT DETECTION 
EXPERIMENTS 



Given the interaction cross section of LKPs with nucleons the expected event rates for 
each experiment can be calculated. A spherical halo model is conventionally used for a con- 
sistent comparison of the different experimental results. Spin-independent (spin-dependent) 
interactions with non-zero momentum transfer to the nucleus demand nuclear (spin) form 
factor corrections in the cross section calculations. Furthermore it is a common procedure 
to normalize the cross sections to the scattering from a single nucleon. In the following we 
describe the details of our event rate calculations for direct detection experiments particu- 
larly with regard to setting limits on the cross sections and couplings. Further information 



can be found in reference 



67] . 



A. Spherical halo model and differential event rates 

The dark matter halo is assumed to be an isothermal and isotropic sphere of an ideal 
WIMP gas obeying a Maxwell-Boltzmann velocity distribution 

f(v,v E )~e , (28) 

where v denotes the velocity of the WIMPs in the rest frame of the earth and ve the velocity 
of the earth with respect to the motionless galactic halo. Vq (= 220 km/ s) is the characteristic 
velocity of the distribution which is assumed to be equal to the galactic rotation velocity 
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v r . There are three contributions to ve- the galactic rotation velocity v r , the velocity of the 
sun with respect to the galactic disc v s and the velocity of the earth around the sun v or b- 
Its main contribution is given by 



VE(t) = V r + V s + V or b COS (3 COS ( 2lT 



t -tp 

T 



(29) 



taking the angle (3 = 59.575° between the earth orbital plane and the galactic plane into 
account. The velocities are given by v r = 220 km/s, v s = 12 km/s and v or b = 29.79 km/s 
respectively, to is the day in a year corresponding to the 2 nd June (tp = 152.5) and T is 
the number of days in a year (T = 365.25). The modulation of the mean velocity is about 
±6.5%. Since we are not investigating time depending properties of event rates such as 
annual modulation effects, we use the mean velocity of the earth (ve) = 232 km/s, averaged 
over one year. The velocity distribution of WIMPs is isotropic in the galactic rest frame. It 
is limited by the escape velocity of the WIMPs (v esc ) from the galactic halo 



\v + v E \ < v e 



(30) 



which yields a maximum WIMP velocity 



V m ax(0,t) = yj V 2 esc - f|(t) (1 - COS 2 6) - V E (t) COS 



(31) 



with 8 being the scattering angle in the galactic rest frame. In our calculations we take v esc = 
544 km/s |64|. We obtain a mean maximum WIMP velocity given by (v max ) = 518.3 km/s 
by averaging over a year and the scattering angle. The maximum recoil energy E^ max of the 
nucleus is given by 



with the kinematic factor 



The differential event rate is 



dR 



Epr ki 



7T 

4 (v E ) 



TP 1 2 

E Rma X = 2 m LKPV max r , 



km LKP m T 
{rrtLKP + m T ) 



erf 



Vmin + (Ve) 

v 



— erf 



v 2 



— e 



(32) 



(33) 



(34) 
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Ro(cr) is the total event rate for a given cross section assuming ve = 0, v esc = oo and 
integrated over recoil energies from Er = to Er = oo. v m i n is defined as the minimum 
velocity leading to a certain recoil energy Er whereas E denotes the energy carried by a 
WIMP with the velocity vq. The ratio ko/ki arises from the normalization of the WIMP 
density distribution. For completeness all relevant formulas are given below: 



2 N, 



o P 



A u m LK p 



av 



E 



R 



E r 



v 



Eq 

h 



-m LKP v 



erf 



7T Vq 



(35) 

(36) 
(37) 
(38) 



Nq is the Avogadro constant and A u is the atomic mass unit. The local dark matter den- 
sity p = 0.3 GeV/cm 3 is taken from 65]. a denotes either the spin-independent or spin- 
dependent cross section which for the case of 71 LKP are given in (fl5|) and (1231) respectively. 



B. Expected event rate: spin-independent interactions 

The spin-independent cross-section given in equation ffl5l) normalized to a single nucleon 
can be written as 

p,n L Pp,n (Q.Q\ 

° scalar ~ ~jjT a scalar ■ {^J) 

Owing to the finite size of the target nucleus, the LKP-nucleus cross section is valid only 
for the case of zero- momentum transfer (q = 0). For non-zero momentum transfer, a form 
factor correction needs to be considered: 

&SI = ^scalar Fgj , (40) 

where F$i is the nuclear form factor. We used Helm's model [7' 

F S i{qr n ) = 3 e 2 , (41) 

qr n 
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with the spherical Bessel function ji, the effective nuclear radius r n = y c 2 + ^n 2 a 2 — 5s 2 
with a = 0.52 fm and c = 1.23v^4 — 0.60 fm, and the nuclear skin thickness s = 1 fm 



67] . 



The spin-independent differential event rate for finite-momentum transfer can then be 
obtained by replacing a with a si in (135]) and thus in (|34"|) . If the WIMP target consists of 
more than one element, the respective abundances of each isotope fi have to be considered. 
The total differential event rate for a specific WIMP target is 



dRsi _ sr^ , dR S I 
dE R ~^ Ti dE R 



The expected number of events for spin-independent interactions for a given experiment can 
then be written as 

N S i= / dE R —fe(E R )MT, (43) 



Qmin 



dE R 



where MT denotes the total detector exposure in kg-days and e(E R ) the WIMP detection 
efficiency as a function of recoil energy. q min and q max denote the lower and upper bound 
of the WIMP-nucleus recoil energy which is considered in the data analysis. We truncated 
the integral at the minimum of the upper analysis limit and the averaged maximum recoil 
energy (E Rmax ). Due to its definition there is a fraction of WIMPs which can give rise to a 
higher maximum recoil energy especially those which hit the detector in a head-on collision. 
However since the differential event rates decrease approximately exponentially the actual 
value of this cut-off only effects the results for WIMPs with low masses (~ lOGeV). We 



used the statistical method proposed in 80J to obtain the XENON10 limits. 



C. Expected event rate: spin-dependent interactions 

Given the spin-dependent cross section ( 12"3"|) and the differential event rate (154]) . the 
expected event rate calculation for axial- vector interactions is similar to the spin-independent 
case. We first normalize the cross section to the scattering from a single nucleon. The nucleon 
spin expectation values are (S p ) = \ and (S n ) = and (S p ) = and (S n ) = |, for a proton 
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and a neutron respectively. Using these values the SD cross section for a single nucleon is 



94 

p,n _ _/^2 2 2 ( AA\ 

^spin — F^p.rflp^n • 1 44 J 



Comparing (|44p to (1231) . the SD cross section with proper normalization can be written as 

p,n _ 3 ^p : n Jn 1 /. p-n 

- 4 "^T J iV + l(5' Pin } 2<7spin ' 1 } 

Similar to the SI case a SD form factor has to be introduced to account for non-zero mo- 
mentum transfer 

ctsd = cr BP mFsD > ( 46 ) 



where Fg D can be written in the form 8J| 



*°°® = Wy (47) 

with the spin structure function S(q). What makes the limit calculations in the SD case 
involved is the fact that S(q) depends on the WIMP-nucleon couplings. For the spin- 
dependent limits for XENON10 we used the form factors obtained using the Bonn A potential 



given m 
evaluated 



In the zero-momentum transfer limit, the spin structure function S(q) can be 



9 7 -1-1 

S (0) = ±^L±±J N (J N + 1)A\ (48) 

7T 



where A is defined as 



Jn 



For finite momentum transfer it is a common procedure to translate the WIMP-proton and 
WIMP-neutron couplings a p and a n into isoscalar and isovector spin couplings and a\ 
using 

a-o — a-p + a n 

a i — Q>p~ a ni (50) 
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so that the spin structure function can be written as 



S(q) = alS 00 (q) + alS u (q) + a aiS i(q) ■ (51) 

Soq, Sn and Sox represent an isoscalar, isovector and interference term respectively. The 
shape of S(q) is determined by the ratio ^ while its magnitude is proportional to a 2 + a 2 . 
Following |78| we consider polar coordinates in the (a p , a n ) subspace: 

a p = a sin 9 

a n = a cos 9 . (52) 

Pure proton and neutron couplings are obtained by setting 9 = 90° and 9 = 0°, respectively. 
Inserting this ansatz into ( 15TT) yields 



S(q) = a 2 ((sm9 + cos9) 2 S 00 (q) + (sm9-cos9) 2 Sn(q)-cos(29)Soi(q)y (53) 



Similar to (143p the number of events is obtained by evaluating the integral 

/•min(q max ,{£; iima:l: )) jni 

Nsd = T,fi dE R ^§°e(E R ) MT . (54) 



Nsd can be re- written in the form 

N SD = A a\ + 2 B a p a n + Ca 2 n , (55) 

with A, B and C being constant for a given WIMP mass. Inserting (j52p yields 

Nsd = a 2 (a sin 2 9 + 2B sin 9 cos 9 + C cos 2 9^j . (56) 

In order to calculate the limits on the WIMP-nucleon SD couplings, for any WIMP mass 
of interest we performed a scan over the angle 9 from 0° — 360°. Thus, as in the spin- 
independent case, a limit on a 2 can be set. Since fl56l) is a quadratic equation the limits are 
expected to be ellipses in the (a p , a n ) subspace. The allowed (a p , a n ) parameter space is 
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restricted to the inner region of these ellipses. 
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